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Abstract 

The Klein-Gordon equation, the MaxweU equation, and the Dirac equation are 
presented as partial difference equations in the eight- dimensional covariant discrete 
phase space. These equations are also furnished as difference- differential equations 
in the arena of discrete phase space and continuous time. The scalar field and 
electromagnetic fields are quantized with commutation relations. The spin-1/2 
field is quantized with anti-commutation relations. Moreover, the total momentum, 
energy and charge of these free relativisitic quantized fields in the discrete phase 
space and continuous time are computed exactly. The results agree completely with 
those computed from the relativisitic fields defined on the space-time continuum. 
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1. Introduction 

In the preceding paper [1], we have presented the Lagrangian formalism for 
the relativistic wave fields in the covariant discrete phase space, as well as in 
the discrete phase space and continuous time. In this paper, we shall choose 
three special examples of relativistic fields [2]. These are the free scalar field, 
electromagnetic fields, and the Dirac spin-1/2 field. Moreover, we shall 
quantize these fields with the usual commutation and anti-commutation 
rules. We generalize variational techniques for the operator- valued second- 
quantized wave fields. The non-singular Green's functions for various differ- 
ence equations are provided [2,3]. Finally, we compute the totally conserved 
four-momentum and charge for three different fields exactly. These compu- 
tations are performed for the wave-equations satisfying difference- differential 
equations in the discrete phase space and continuous time only. (We do not 
consider wave fields in the covariant discrete phase space to calculate the 
totally conserved quantities for some physical reasons.) The results of these 
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computations are identical to those calculated from the usual relativistic 
quantum theory of the free fields in the Minkowski space-time. One may 
wonder about the utility of such a complicated, alternate formulation of the 
quantum theory of free fields! Actually, the present treatment of the free 
quantized fields in the discrete phase space and continuous time is just a 
prelude to the more exciting investigations of the interacting fields and the 
^-matrix in the following paper III. 



2. Definitions and notations 

We use physical units such that h = c = i = 1. (Here, ^ is a characteristic 
length.) All physical quantities are expressed as dimensionless numbers. 
Greek indices take from {1,2,3,4} and roman indices take from {1,2,3}. 
The summation convention is followed. The Minkowski metric is denoted 
by Vnf the signature of the metric is +2. We denote the set of all 
real numbers by M and all non-negative integers by N. A bold roman letter 
indicates a three-dimensional vector. The equations in the covariant discrete 
phase space are denoted by (■■A), whereas the equations in the discrete phase 
space and continuous time are labelled by {--B). 

An integer n'^ is associated with a phase space circle of radius ^/2nM^\^ 
for each G {1,2,3,4}. Let a function be defined by / : — >^ R (or 
C). The right-difference, the left-difference, and the weighted mean 
difference are defined respectively by: 

A^/(n):=/(..,n'^ + l,..)-/(..,n^..), (li) 
A;/(n):=/(..,n^..)-/(..,n'^-l,..), (lii) 

A#/(n) := (1/V2) [V^^^/(..,n'^ + 1, ..) - f{.., n'' - 1, ..)] . (liii) 

Let us work out an example to elaborate the definition (liii). Suppose that 
ki is a real number and H^i (ki) is a Hermite polynomial. (See Appendix L) 
In that case 



Afeni(fci) := A* 



(z)^^e-(fei)'/^g„i(fci) 
7rV42"V2y(%! 

H^i+iiki) + 2niHni-i{ki] 



.(nl+l)g-(fci)V2 



27rV42"V2^(7^! 

i)"^e-('=i)'/2iJ^i(/t^ 



(2) 



(ik. 



7rV42nV2y(^! 



(ifei)^„i(fci 



In the case of 



n 



(/)(n^, n?, n^, n^) 



is an operator-valued function 



over a domain of N , we adopt the same definitions as in (li,ii,iii). 
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3. The second quantization of a free non-hermitian scalar field 



Let p, Pn be five linear operators over a non-separable Hilbert space [4]. Let 
denote the corresponding adjoint or hermitian-conjugate operators. 
The linear operator- valued Lagrangian function L of ten operators is defined 
to be [5]: 

L{p, P^\ Pix, pI) ■■= -{n^^'plpu + mW) , 



dL{..) 
dp 


2 t ^-^(-O 2 
-pp\ = pp, 






dL{..) 
dpfi 


nv \ dL{..) _ 
'/ Pi/ 1 _ t ~ ' P'^ '> 

dpi 




(3) 


d^L{..) _ 
dpdp^ 


d^L{..) , d^L{..) 

dp^dp ^ ' {dpf 






d^L{..) 
dpidpu 


dpudpl ' ' ^p^^p^, 


= ^'^("^ - 

dpidpt 





Here, p > is the mass parameter, "I" is the identity operator, and "O" is 
the zero operator. The linear operators do not commute necessarily. There- 
fore, the order in which they appear must be preserved. The j -th partial 
derivatives of L in (3) with respect to p, p^P/j,, p'^^^ are all zero operators for 
j > 3. Assuming that the Euler-Lagrange equations (38A,B) of paper I are 
valid for the operator- valued functions p = (f){n) and p = (j){n,t), we obtain 
[2] from (3) 

77'^^A*A*<^H-MV(n) = 0, (4A) 
<5«^A#A#<^(n, t) - (dtfHn, t) - ^^(n, t)=0. (4B) 

(Strictly speaking, is a section of the fibre bundle [6] of linear operators 
over the base "manifold" N"^ or x M. ) 

Now we shall derive the operator difference conservation equations. Adapt- 
ing the equations (A.II.2), (A.IL3), and (A.II.4A) of paper I to the 
Lagrangian in (3) and noting that operator ^g^g'J etc. are proportional 
to the identity, we obtain relativistic conservation equations: 





r9L(..)i 




. dpt, _ 



dL{..) 



+- 



dpWp dpdp^ 



-0t(n).A*0-[A*0t].<^(„)} 
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1 

+ 2 



a2L(..) ^ d'^L{.. 



dpadpt, d Pad pi 



A 



+ 



AJ(1) 



aL(..) 



0(n) 



dp^dp dpdp^ 



dp\.. 

dPul dpi 
</.t(n)0(rz) 



2 [Opt^p,, d Pad pi 



0. 



(5A) 



Here, (h.c.) stands for the hermitian-conjugation of the preceding terms. 

Since the Lagrangian function L is a second degree polynomial, the 
relativistic conservation equation (5A) is exact and no additional terms in- 
dicated by . . . are necessary. 

Moreover, the last curly bracket in (5A) is exactly the zero operator. 
Using (5v) of paper I and (6i), after a long calculation, the relativisitic 
equation (5A) yields the exact conservation equation: 





\ dL{..) 1 




dpu\{..,n- -!,..)_ 



A*</> + ^-[A#</,]|(.._„._,_..) + (h.c.) 



(^P^-.n-l,..) 



•A^ 



(6A) 



-p\<t^\..,n- - 1, ..) • 4>{n) + (t>\n) ■ ct>{..,n- - 1, ..)) 



0. 



In a similar fashion, we can derive exact difference-differential equations 
(which are equivalent to the relativistic equations): 



A;,r„^ + 9ir„^ = o, 

A6T| + 5tri = 0, 



(6Bi) 
(6Bh) 
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dL{..) 



dpb\.. 



(A,V)|(..,..-i,..)-A|</. + A*</,t-(A, 



l'^)|(..,n''-l,.) 



+ 



idt^^)\{..,n''-l,..) ■ dt(t> + 9t<^^ • (5t<^)|(..,„6_i,..) 



2 Ut (.., - 1, ..) . </,(n) + </.t (n) • - 1, ..)] , 



(6Biii) 



5P4|.. " 



V 2 [ dpa\.. 
Ti{n,t) := • dtcP + (h.c.) - L(..) 



(6Biv) 

(6Bv) 
(6Bvi) 



It is instructive to compare and contrast the equations (39Bi,ii,iii)iv,v,vi) 
of paper I with the equations (6Bi,ii,iii,iv,v,vi). The first set of T(f{Ti,t) 
satisfy possibly approximate conservation equations whereas the second set 
satisfy exact equations. Neither of these Tl}{n,t) obey exactly the tensorial 
transformation rules (34B). However, the relativistic total four-momentum 
components P^'s and the invariant total charge Q can be elicited from 



oo(3) 



n=0 



0L{..) 
dp4\.. 



• A*</> + A*<; 



dLj..) 
dpi 



4|. 



>\t=0 



(7Bi) 



H = - 



oo(3) 



( ^P^\- dpl^. 



-L{-)\. 



'\t=o 



(7Bii) 



Q = -^e^lf^. <P{n, t) - </>t („, t) ■ ^ 



dp. 



4|- )\t=0 



(7Biii) 



Note that above relativistic equations are exac-t and no additional terms de- 
noted by . . . are necessary. We furnish a general class of exact solutions [2] 
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("the plane wave superposition") of the (generalized) Klein-Gordon equa- 
tions (4A,B) in the following: 



(j){n) 



[2a;(k)]-V2 a(k) 



+ 6t(k) 



,1=1 



d^k 



=: (l>-{n,t) + cj)+in,t), 



n ^n.(%) 











3 




|a(k) 













(8A) 



(SB) 



(k) := (fci, k2, ks), d^k := dki,dk2,dks : 



-k4 = u{k) := +^(kiY + {k2f + {k^Y + lJ? > , 

[-(A;^)V2]-i^„M(fc^) 

UmIM-W (^)l/4 2nM/2^(7^ • 



(9i) 

(9ii) 
(9m) 



Here, the indices jj, and j are noi summed. Moreover, HnP-ikfj) stands 
for a Hermite polynomial. (For the properties of orthonormal complex- 
polynomials (,ni^{kfji), see equation (2) and Appendix I.) The functions "a" 
and "6" are some sections of the fibre-bundle [6] of the linear operators 
over the base manifold (the momentum-space). The operator- valued 
improper integrals 

j at(k)a(k)d^k, j a(k)ot(k)d^k, j 6t(k)6(k)d3k, j h{\i)h\\<i)d^\L, etc. 



should converge in certain sense for the existence of (8A,B). There are more 

restrictions on these operators which follow from the quantum theory. These 
are the following canonical quantum rules to be imposed on the operators 
o, a^, 6, 6^ : 
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[A, B] := AB-BA, 

3 

d^{k - k) := 6{ki - ki) 5{k2 - fe) 5{k^ - aJs) = n '^(^i - ^i) ' 

j=i (10) 

[a(k), at(k)] = [6(k), 6t(k)] = <53(k - k) I (k) , 

[a(k),a(k)] = [at(k),at(k)] = [6(k),6(k)] = [6t(k), 6t(k)] =0. 

Here, 5{kj — kj) denotes a Dirac-delta distribution function, I (k) stands for 
the identity operator, and "0" denotes the zero operator. The Unear opera- 
tors a(k),a^(k) are called the destruction and creation operators for parti- 
cles (or field quantas). The particle and anti-particle vacuum is denoted by 
the Hilbcrt vector |V'o)- The particle and anti-particle (occupation) number 
operators (k) , (k) are defined by and satisfy the following equations: 

N+(k) := at(k)a(k), Ar-(k) := 6t(k)6(k) , 

a(k)IV'o) = WIV-o) = iV+(k)|V'o) = iV-(k)IV'o) = |0) , 

(V'o|^o) = l, (0|0)=0, 

[iV+(k),a(k)] = -53(k-k)a(k), 
[Ar+(k),at(k)] =5'^{k-k)a^k), 

(11) 

[iV-(k),6(k)] = -(53(k-k)6(k), 
[Ar-(k),6t(k)]=<53(k_k)6t(k), 

[Ar+(k),6(k)] = [iV+(k),6t(k)] = 0, 

[iV-(k),a(k)] = [iV-(k),at(k)] = 0, 
[iV+(k),iV-(k)] = 0. 

The eigenvalues of Ar+(k) and A'^~(k), the so called occupation numbers, 

take values from N := {0,1,2,3,...}. Therefore, the particles and anti- 
particles of the quantized scalar field obey the Bosc-Einstein statistics. 

The covariant commutation relations which follow from (8A) and (10) 
are the following (see Appendix II): 

[<p-{n),<p-{n)] = [<^+(n),0+(n)] = [(0- (n))t, (<^-(n))t] 

= [(</.+ (n))t,(0+(n))t]^O, 
[0-(n),(<^-(n))t] = (i/27r)A+(n,n;/x)I, 

12A 

[,^+(n),(<^+(n))t] = (V27r)A_(n,n;/x)I, 
[^(n),</.(n)] = [(0(n))t,(<^(n))t] = O, 
[</-(«), (0(n))t] = (i/27r)A(n,n;M) I. 
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Moreover, the covariant commutation relations in the difFerence-difFerential 
representation are: 

[nn,t),r{n,t)] = [<^+(n,t),0+(n,?)] = [(<^-(n,t))t, (^-(n,F))t] 

= [(<^+(n,t))t,(</,+ (n,r))t]^0, 
[0-(n,t),(0-(n,?))t] =^A+(n,t; n,t- 

[0+(n,t),(,/>+(fi,t))t] = iA_(n,t; n,f; 
[(/>(n,t),</.(n,r)] = [0(n,t)t,(0(fi,r))t] = o, 

(12B) 

[</.(n,t),(0(n,t))t] = iA(n,t; n,t; 
[</.(n,t),(</.(n,i))t]|^-^^ = for n^^n, 

[dt<P{n,t),i<P{n,t))%=t = -i {\{^n^n^^=- -^nnl, 

[dt<l){n,t),{di^{iy,t))\=t = ^ for "T^fi- 

The last three commutators in (12B) resemble the three fundamental pos- 
tulates of quantum mechanics, namely 

[Q„, Qb] = 0, [Pa, Qb] = -iSab I, [Pa, A] = . 

Now we shall compute the total three-momentum components Pj , the to- 
tal energy H, and the total charge Q from the equations (7Bi,ii,iii). (See the 
comments at the end of Section 5 of paper I for not considering equations la- 
belled A.) Elaborate computations are explicitly performed in Appendix III. 
We summarize the results in the following equations: 

Pj = J [N+{k) + N-{k)] kjd^k , (13i) 
H = -Pi = y^[iV+(k) +iV-(k) + (53(0)1 (k)]w(k)ci3k , ((ISii) 
Q = e j [iV+(k) - Ar-(k)] d^k. (ISiii) 

9? 

These results are identical to those derived from the usual relativistic quan- 
tum theory of a free non-hermitian scalar field [7] in the (flat) space-time 
continuum. The equation (13ii) shows that the divergence of the null-point 
energy cannot be remedied by the discrete phase space approach. 

We now use the commutators in the equation (10), (11), the field in 
(8B), and the conserved operators in (13i,ii,iii) to derive: 

[Pj,^{n,t)] = -iA*<f>{n,t) , 

[H,(t){n,t)] = -idt(t>in,t) , (14) 
[Q,cf>{n,t)] = -ecf>{n,t). 



8 



The above equations prove that the operators Pj,H, Q are the generators for 
the space translation, the time translation, and the gauge transformation. 
(We have put here (^3(0)0 = 0.) 

Conservations of the total momentum components Pj and the total 
charge Q can be proved alternatively by the commutation relations: 

' (15) 
[H,Q] = 0. 



4. Quantization of free electro- magnetic field 

The Lagrangian function L for the electro-magnetic field is chosen to be 
the following second degree polynomial function of the twenty [5] linear, 
self-adjoint operators and Uni,: 

= -{l/2)[S^f'S'^yacybd - S^'^yoAyM + VAaVAb) + (y44)'] , 



dyti ' dypr ^ ' dyab 

dL{..) dL{..) dL{..) 

= Via , = UoA , = —2/44 

oy-ia oyaA oy^A 



(16) 



{dyf,u)idypr) 
d'^L{..) d'^L{..) d'^L{..) d'^L{..] 



{dVabY (dyoAf {dyiaY {dy^A)'^ 

d^L{..) _ d^L{..) 
{dya4){dycd) {dy4a){dycd] 



= -I. 



= 0. 



The third and higher order partial derivatives of L are obviously all zero 
operators. 

The Eulcr-Lagrangc equations (36A,B) of paper I extracted from (16) 
with y^ = Afj,(n) etc. yield: 

ni^-A*A*A,{n) = 0, (17A) 
S-b^*A*A„{n, t) - {dtfA^in, t) = . (17B) 

These equations are further augmented by the Lorentz-gauge constraint 
on the allowable state vectors l'^) (in a "Hilbert space" with indefinite met- 
ric): 
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(V|A*^'^(n)|^)=0, 
{'iP\A*A\n,t)+dtA\n,tM)=0. 



(18A) 
(18B) 



The Maxwell's equations (17A,B) and the Lorentz-gauge constraint (18A,B) 
are preserved by the restricted gauge transformations involving a hermitian 
operator J7: 



Aj{n,t) = Aj{n,t) - A*n{n,t) , 
A4n,t) = A4{n,t) - dtn{n,t) , 

5'^'AtA*^{n,t)-{dt)^n{n,t) = 0, 

{ilj\A*Ai'{nm=0, 
{i^\A*A''{n,t)+dtA\n,tm=0. 



(19A) 
(19Bi) 
(19Bii) 
(20A) 
(20B) 
(20A) 
(20B) 



By the Lagrangian (16) and the equations (39Biv) and (39Bvi) of paper I, 
we obtain 

oo(3) 



n=0 



|t=0 



oo(3) 



H = -P, = (1/2) [5^\AtAn-{AtA,) + idtA'^).{dtA, 



n=0 



|t=0 



(21Bi) 
(21Bii) 



The above relativistic equations are exact and no additional terms are nec- 
essary. The "plane wave" decomposition of the four-potential operator 
An{n,t) is given by [2]: 



A^{n,t) = 4(n,0 = I [2Kk)]-vJa^(k) 



-iut 



n ^n.(fc.) 



e-* d3k=:A-(n,t) + A+(n,t)^ 



A*A^(n,t)|t=o = i I [2z.(k)]V2fcJa^(k) 



j 



-at 



<(k) 
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H ini {kj 



-«i(k) 



(22B) 



The canonical quantization rules are assumed to be: 

K(k),a,(k)] = [at(k),at(k)]^0, 
K(k),at(k)]=r?^,<53(k_k)I(k). 



(23B) 



The covariant commutation rules, which follow from (22B), (23B), and 
(A.II.5B), are summarized below: 

[A+{n,t),A+{n,t)] = [A~{n,t),A-{n,t)] = 0, 

= iVnuD+{n,t; n,t)I, 
[A^i^,t),A-{n,t)] = iri^,i,D_{n,t; n,t)I, 

[A^{n,t),A^{n,t)]=ir]^^D{n,t; n,t)l, (24B) 
[A^{n,t),A^{n,t)]^i^^ = n/fi, 

[dtA^{n,t),A,{n,t)]ii=t = -iVnuS^-I, 
[dtA^{n,t),dfAi,{n,i)]^f^^ = for nT^n. 

The computations of the total momentum-energy from (21Bi,ii) and (22B) 
yield: 

Pj = (1/2) J ry'^'^[aj,(k)a,(k) + a^(k)at (k)] kjd^k , 

r (25B) 
H = -P^ = (1/2) J <^[aj,(k)a,(k) +a^(k)at(k)]i^(k)d3k. 



We can choose a special gauge so that (25B) simplifies considerably. We 
assume the condition 



(V^lfcXWIV') =0. 



(26) 



This is a sufficient condition for the satisfaction of the Lorentz-gauge condi- 
tions (18A,B). We introduce a special restricted gauge condition (see (22B)) 
with the help of the hermitian operator-valued function 



-at^(k) 



-iut 



(27B) 



d-^k. 
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Under this gauge transformation, the new field operators aju(k) (correspond- 
ing to the field A^(n, t)) undergo the following transformations: 

a^(k) = G^(k) - [,y{k)]-^k4kf,a\k)] , (28i) 

(V|fc^a^(k)|V')=0, (28ii) 

(V'|a4(k)|V') =0. (28iii) 

Thus, by (28iii) the temporal component 04 (k) drops off. Next, let us 
consider the "orthonormal" tetrad [8] e|^^^(k) (also see equation (23) of 
paper I) which in general satisfy: 

r?(^-)ef,)(k)e^^)(k) = <-, 

(29) 

«(A)(k) := a^(k)e{\)(k) , 
a'^(k) = ef,)(k)aW(k). 

We can choose prudently (for z/(k) > 0) two of the tetrad vectors by the 
following: 

(e;3)(k), e23)(k), e3 (k), (k)) := [u{k)]-Hh,k2,h,0) , 
<4)(k):='^f4)- 

The choice of the other two vectors e^-^j(k) and e^2)('') arbitrary up to a 
two-dimensional orthogonal transformation. 
The condition (28ii), by (30) yields 

fc^ef,)(k)(V'|aW(k)|^) = k,e^^^{k){^\a(^){km 

= Kk)](^|a(3)(k)|V')=0. 

Thus, the expectation value of the longitudinal component {tp\a^^\k)\ip) 
vanishes. Dropping circumflexes in the sequence, we obtain from (25B), 
(28iii), (31), and (29), the simplified versions of the expectation values of 
the total momentum-energy as: 



{m\^) = T. /(V'|iV(A)(k)|V)M'k, 

^=1 R3 

(mm = E / (V'iA^(A)(k) + {i/2)s^{o)i^x){k\^)u{k)d'k. 



A=l 



Only the two degrees of (linear) polarization (A G {1,2}) contribute to the 
total momentum-energy of the photons. The equations in (32B) are identical 
to those obtained by the usual relativistic theory [7] . 
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5. Quantization of free spin- 1/2 field 

The 4x4 Dirac matrices satisfy [5]: 

, (33) 

The Dirac bispinor field p = ^(n) or p = ^p{n,t) is a 4 x 1 column vector 
of operators in the second quantized theory. We denote 

• t 4 

t (34) 
[pp] = PP- 

The Lagrangian function for a massive, spin-1/2 field operator p is taken 
to be: 

Hp, p; Pu-, h) ■= -(.^/^){pi^Pu- - PhI^p) - mm m>o, 

^I^ = {l/2)p,r-m~p, ^ = -(1/2)7 V^-mp, 
dpdp dpdp ' dpdpij, dp^dp 

^ = ^ = (l/2)7^ 

9p^9p dpdp^ 

d'^L{..) d'^L{..) d'^L{..) d'^L{..) 



(35) 



dp^dp^ dp^dpy dpij,dp„ dp^dpy 



= 0. 



The triple and higher partial derivatives of L are all identically zero opera- 
tors. 

The Euler-Lagrange operator equations from (35) and equations (38A,A), 
and (38B,B) of paper I are furnished by: 

7''A*V(n) + mV'(n) = , (36A) 

[A*V'(n)] 7^" - rwipin) = , (36A) 

7^' A|V(n, t) + 7^atV(n, t) + mV'(n, t) = , (36B) 

[A*V'(n, t)] 7^' + [5t^(n, t)] 7^ - "iV'(n, i) = • (36B) 

By the equation (35), and equations (A.II.5A) and (A.II.6Bi,ii) of pa- 
per I, we derive: 
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At 



dL{..: 



dL{..) 
dpu 



A*A*V + (h.c.)-A*L(..) 



-Vi(n) • A* 



[dpdp, 



dpdpj, 



A*V 



+ A* 



A*^| 



dpudp 



dp„dp 

+0=: A,r;(n) = o, 



A^V' - (A*A*v;)- 



dpudp 



A,T^{n,t) + dtT^{n,t) = 0, 
Af,T!{n,t)+dtTi{n,t)=0, 



dLj..) 

dp^i.^nb-l,..) 



Af^ 



+?^.(A#^)|(..,,._i,..) + (h.c.) 



dpb\.. 



-s: 



1 .r. 



n" - 1, ..) • j'^A*^ + ^(n, • nAti^)i(..,n^-i,..) 

+^(.., - 1, + Vi(n, 07^(5tV)|(..,n<.-i,..) 

-(A*^)|(..,n''-i,..)-7^V'(n,t) 

-(A*V;)7^^(-, - 1, ..) - (5tVi)|(..,n<'-i,..) • ■ V'Cn, 
-(9tV')-7'-V'(",n''-l,..) 

+m(V!(.., - 1, ..) • ip{n, t) + i^{n, t) ■ - 1, ..) 



14 



dp4\. 



dp4\. 



(37Biv) 



T!{n,t) :-- 



dL{..) 



%|(..,n'-l,..) 



(37Bv) 



ri(n,t) :-- 



5p4|.. dp4.\.. 



L{..\. 



(37Bvi) 



The above equations are exact. Note that the tricky ordering of the operators 
in the above equations is crucial. 

The relativistic total momentum-energy and the charge are given by (see 
equations (7Bi,ii,iii) and (37Biv,vi)): 



00(3) 



Pj = m Y: [A|V^-V(n,i)-V^(n,t).A*V]|t=o, 



(38Bi) 



n=0 



00(3) 



n=0 



H = -Pi = {i/2) [V-^Cn, t) ■ dti^ - ■ V'(n, t)]|,=o ' (38Bii) 

(38Biii) 



00(3) 

Q = eY V'Hn,0)-V(n,0) 

n=0 



A class of "plane wave" solutions of the Dirac equations (36A) and (36B) 
is provided by (see reference [2] and Appendix I): 



+/3t(p)v,(p)(^n^^^) 

^(n)= /[m/E(p)]V2/^ 

R3 1^ r=l 

=: ip~^{n) + ip~{n) , 



(39Ai) 



aj(p)ur(p) n^'»''(^^/^) 

Up 



(39Aii) 
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^in,t) := |[m/^(p)]V2|^ 
=: V-(n,i) + V+(n,i), 



a. 



-iEt 



(39Bi) 



ras ^ '■=1 



=: ^^+(n,i)+V;-(n,i), 



-iEt 



at(p)u,(p)(nen.fe))e^''* 

3 (39Bii) 



xi-'p 



P = -Pi = E{p) := +\l5"-^PaPh + m2 > , 
ut(p)us(p) = vj(p)vs(p) = [E{p)/m] 5rs , 
u]:(p)vs(-p) = vj(p)us(-p) = 0. 



(39Bm) 



The canonical quantization rules for a spin- 1/2 field operator are fur- 
nished by the anti-commutators: 



[ar{p),as(.p)]+ 
[ar{p),Ps{p)]+ 
K(p),aJ(p)]+ 



AB + BA, 

[Pr{p),Ps{p)] + 

[at(p),«t(p)]+ = [/jt(p),^t(p)]^^0, 
[aUp),Plip)U 

K(p),/3l(p)]+ = [alip),Psm+ = 0, 

[Mp),Plm+ = SrsSHp-p)^- 



(40) 



The particle and anti-particle occupation number operators are defined 



by: 



iV-(p) 



4(p)"r(p) , 
PKp)Prip). 



(41) 



Here, the subscript "r" is not summed. The occupation number operators 
N~{p) and iV+(p) take eigenvalues from {0,1}. Therefore, this quantiza- 
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tion is compatible with the Fermi-Dirac statistics. The particle and anti- 
particle vacuum state I'^o) is characterized by: 

MpMo) = MpMo) = |0), 

(42) 

(V'o|V'o) = l, (0|0)=0. 

We can derive the covariant quantization rules for a spin-1/2 particle 
field by equations (40), (39A), and (39Bi,ii). These are provided by the 



following anti-commutation relations: 

[V'-(n),Vi-(n)]+ = [V'+(n),Vi+(n)]+ = 0, (43Ai 

[V-(n),V'+(n)]+ = [^p-{n),^p+{n)]+ = 0, (43Aii 

[V'(n),^(n)]+ = [i/i(n),Vi(n)]+ = 0, (43Aiii 

[i;-{n),i>^{n)]^ = iS+{n,n)l, (43Aiv 

['0+(n),'0~(n)]+ = -iS-{n,n)l, (43Av 

[ip{n), i^{n)]_^ = —iS{n, n) I , (43Avi 

[^-(n,0,i/i-(n,O]+ = [V'+(n,t),^+(n,i)]+ = 0, (43Bi 

[V'-(n,t),V'+(n,?)]+ = [V;-(n,i),V^+(n,F)]+ = 0, (43Bii 

[^(n,t), V(n,f)]+ = [i/i(n,t),^(n,?)]+ = 0, (43Biii: 

[V'-(n,t),^+(n,f)]^ = -i5+(n,i;n,?)I, (43Biv 

[V'+(n,i),^-(n,?)]+ = -i5_(..)(n,i;n,?)I, (43Bv 

['i/'(n,t),Vi(n,f)]^ = -iS'(n,i;n,?)I, (43Bvi 

[V'(n,0),V^(n,0)], =(53-1. (43Bvii 



Here, the Green's functions S (n, n) and S (n, n, f ) are borrowed from 
Appendix II. ^""^ ^"^ 

Now, we shall compTitc the total momentum, energy, and charge by the 
equations (39Bi), (38Bi,ii,iii), (40), and (41). After laborious calculations, 
similar to those in Appendix III, we obtain rather neat results: 

P3 = jl J[N-{p) + N+{p)]pjd'p, (44Bi) 
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H=-P4 = i2 I i^riP) + ^r+(P) - '^'(0) Ir(p)] E{p)d'p , (44Bii) 

Q = e^ I [iV-(p) - iV+(p) + <5=^(0) I,(p)] d^p , (44Biii) 

In the case of the electron-position field, we choose A^^(p) as the number 
operator for electrons, and the charge parameter e = — \/47r/137 ; in contrast 
to the charged scalar particles in (13iii). The results (44Bi,ii,iii) coincide 
exactly with the usual relativistic quantum theory of a free spin- 1/2 field in 
the (flat) space-time continuum [7]. 



Appendix I: Hermite and related complex polynomials 

The definitions and the useful formulae for the Hermite polynomials [9] are 
provided here (fc G M; n G {0, 1, 2, . . .}) : 

iI„(fc):=(-l)V^^(e-'=^), (A.I.1) 

Hni-k) := (-Ifi^nW , (A.I.2) 

d^Hnik) ^, dHr,{k) 



- 2k — ^ + 2nHn{k) = , (A.I.3) 
dk 

2nHn-i{k), n>l, (A.I.4) 



dk"^ dk 
dHnik) 



dk 

Hn+i{k) = 2kHn{k) - 2nHn-i{k), n>l. (A.L5) 

The complex, orthonormal polynomials £,n{k) satisfy the following definition 
and equations: 

1^1/4 2«/2^ -^(^)"^--(^)' C-n-iW:=0, fornGN, (AL6) 



U-k) = i-lTUk) = Uk) , (A.I.7) 

d^Uk) ^ ^_^2 ^ 2n + l)Cn{k) = , (A.I.8) 

+kUk)=iV2^U-iik), (A.I.9) 

A*Uk) = ikUk) , (A.I.IO) 



d^n{k) 
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(-l)"e 



7rV4 2"V2n! 

n 



(2fc)2" + ^(-iy22'^-^(2j - 1)!! fe2(n- 



-i) 



7rV4 2"+VV(2n+l)! 
(2fc)2n+i +^^^(_l)j22"+i-J(2j - 1)!! (^^"^+ ^) 

((2i-l)!! := (2i-l)(2i-3)...5.3.1 .) 



2i 



Co(A;) = 



7rV4 ' 



TT 



1/4 



6n(0) 



(2n-l)!! ^{2ny. 



7rV4 7(2^ 7rV4 2"n! ' 

C2n+l(0) = 0, 



exp {[(t^ + A;2)/2] + i^tk} = V 



2"/2 exp { - [(fc - v)lt}U{k + P)/ V2 ] 



.V4^ 



n 



^ ^ ^ , 

ik-k)Y, Uk) Uk) = i^iN + l)/2 



n=0 



CXJ cxj 

j im{k:)^nik)dk= j fm{k)fn{k)dk = 6mn, 



Y.Uk)Up) = 5{k-p). 

n=0 
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oo oo 
n=0 n=0 

oo 

C„(A;) = -^ / U{x)e'^-dx. (A.L21) 

— oo 



Appendix II: Non-singular Green's functions 

The relativistic invariant Green's functions [2] for the finite difference and 
the difference-differential Klein-Gordon equations (4A) and (4B) respec- 
tively are given by (see Fig. 1): 




Cf 



FIG. 1 The complex A;^ -plane. 
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(Note that in our signature k'^ = — /C4.) The Green's functions involving 
the closed contours in Fig. 1 are called homogeneous, whereas the Green's 
functions involving the open contours are called inhomogeneous. Assuming 
appropriate uniform convergences of the improper integrals in (A.II.1A,B) 
and using the equation (A.I. 18), we derive that 



^M-A# A# A(„) (n; n; fi) - //^ A(„) (n; n; m) 



/ / 

Jm3 Jc{a) 



for the inhomogeneous, 
for the homogeneous; 



(A.II.2A) 



^^■^ A* Af A(,) (n, t; n, t; fi) - 0^ A(„) (n, t; n, t; //) 
-1 



= -(27r)-^ - 



ex.p{ik4{t — t))dk'^ 



C{a) 

^ ^ f~'^nn^(^~0 for the inhomogeneous. 







for the homogeneous; 



(A.II.2B) 



Consider the particular homogeneous Green's function A(n, t;n, t; ^u) in 
(A. II. IB) which is associated with the closed contour C in Fig. 1. Performing 
the closed contour integration in the complex fc*^ -plane, we obtain 



A{n,t;n,t;n) = 

A{n, t; n, t; iJ,) = for 117^11, 

[dtA{n,t;n,t;fi)]^i^^ = -6l^, 
[dtdt A{n, t; n, t; fi)]^f^^ = . 



[u;"^sin(u;(i-f))]d^k. 



(A.II.3B) 



There exist the following relations among the Green's functions: 
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A(n, t ; n, f; /x) = A+ (n, t ; n, F; //) + A_ (n, t ; n, t; jj) , 



A_(n, t; n, t; fi) = A+(n, t; n, t; fi) , 
A(n, ?; n, t; fi) = - A(n, t; n, ?; ^) , 

AF(n,i;n,?;//) := -6l(t - f)A+(n, n, f; ^) (A.II.4B) 
+6l(f-t)A_(n,t;n,f;//), 

AF(n, t; n, i; //) = Apin, t; n, i; /x) , 
0(i):=Q)(l+t/|i|) for tT^O. 

We define the Green's functions for the massless case /i = by: 

(n; n) := A(„) (n; n; 0) , (A.II.5A) 

(n, n, i) := A(„) (n, i; n, f; 0) . (A.II.5B) 

In the second quantization of the spin- 1/2 fields, we encounter the 
4x4 matrix- valued Green's functions [2] S^^) ip-'i ^5 S^^) (n, t\ n, t; m). 

These are defined by: 

S{a) {n; n;m):= {j^'A* - m I) A(„) (n; n) 

irj"^PaP/3 + m'^)-\ij>'p^ -ml) 



(a) 



t/=l 



(A.II.6A) 



S{a)i^,t;n,t;m) := (7-'A| + 7^5* - mI)A(„)(n,t; n,?;m) 



= (27r)-^ J i J {t^PaPp + m')-\i^^p^ + i7>4 - ml) 

6=1 J J 



(A.II.6B) 



Here, m > is the mass parameter, 7^ are Dirac matrices, and C(a) are 
contours in the complex -plane (exactly similar to those in Fig. 1). 
We can prove that 

{—6t~ I for the inhomogeneous, 
(A.II.7A) 
for the homogeneous. 



22 



(y A* + + m l)S^a) (n, t; n, t; m) 

{—S'^^6{t — t)l for the inhomogeneous, (A.II.7B) 
for the homogeneous. 

There exist Unear relationships 

S{n; n; m) = >S'+(n; ra; m) + S-{n\ n; m) , (A.II.8A) 

^(n, t; n, t; m) = S+{n, t; n, t; m) + ^-(n, t\ n, t; m) , (A.II.8B) 

SF{n,t;n,t;m) := —9{t — T)S+{n,t;n,t;'m) 



+9{t-t)S-{n,t;n,t 



(A.II.8C) 



;m 



Appendix III: Total momentum, energy, and charge of the scalar 
field 

Wc shall compute here the total momentum Pj , the total energy H = —P4, 
and the total charge Q from the non-hermitian scalar field operator ^(n, t) 
in the equations (7Bi,ii,iii). We have from (SB), 



4>{n,t) = |[2a;(k)]-V2|a(k) 

]R3 

+6t(k 

|[2a;(k)]-V2Jat(k) 



L j=l 



—iu!t 



L j=l 



+6(k) 



-fet(k) 



>d-^k, 



dtcP{n,t)\t=o = -i / a;(k)[2a;(k)]-i/2<ja(k 



-5t(k)nc^(fc£) 



(A.III.l) 
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We shall first compute the total charge Q from (7Biii) and (A.III.l) for 
the sake of simplicity. 

oo(3) 

Q = ieY^ [0t(n,t)a,0(n,t)]|,^o + (h.c.) 

n=0 

oo(3) 



= («/2) E / / i^/^) 

n=0 „Q „!! 



1/2 



+ 



^2 1 n n ^n4„^ - a^st n n 



+ (h.c.) . 



Carrying out the triple sum J2 with help of the completeness relations 
(A.I.19) and (A.I.20), we obtahi 

Q = (e/2) y" J{Q/uy/'^[a''a-bb^)S^k-k) 

+{ba - a^b^)S^{k + k)] (fkd^k + (h.c.) 
= (e/2) I {[at(k)a(k)-6(k)6t(k)] 

R3 

+ [5(k)a(-k) -at(k)6t(-k)]}(i3k+ (h.c.) 
; j [at(k)a(k) - b{k)b\k)] d^k . 



(A.III.2) 



Using the definitions (11) and the commutators (10), we get 
Q = e l[N+{k)-N-{k)-6^{0)l{k)]d^k. 



The last divergent term in the above equation could have been avoided by 
modifying [10] the Lagrangian (3). Thus we derive the equation (13iii). 

Now, we shall compute the total momentum components Pj from (A.III.l), 
(7Bi), and (10). The result is 

oo(3) 



n=0 

= (1/2) E / y kjiO/u)'/' 



n=0 , 



\ e j / \ e j 
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d^M^k + (h.c.) 



= (1/2) J fc,{[at(k)a(k) -5(k)6t(k)] - [at(-k)6t(k 

+6(-k)a(k)]}d3k + (h.c.) 

A;j[at(k)a(k) + 6^(k)6(k) + 6^{0) I(k)]d^k 



(A.III.3) 



Here, we put S^{0) J kj !{]<.) d^]s. = in the sense of the Cauchy-Principal- 

Value and thus derive the equation (13i). 

Finally, we calculate the total energy H by the equations (A.III.l), 
(7Bii), (10), and (11). We obtain 

H = -Pi 



oo(3) 



= [(^'^''(A*0t.A*</,) + 5t<^t.a^<^ + ^2^t(n,i)<^(n,i) 



n=0 



|t=0 



(1/2) E 



-1/2 



(ata + 66t)53(k_k) 

- {S^^kak^ + UU- /x2)(at6t + ha) (^^(k + k)l d^kd^k 



(A.III.4) 



0J[ 



(k) fat(k)a(k) + 6(k)6t(k)l d^k 



By the above equation, the commutation relation (10), and the number 
operators in (11), the total energy in (13ii) is obtained. 
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FIG. 1 The complex /c^ -plane. 



